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Abstract 

Liquids in contact with solids are submitted to intermolecular forces making liquids 
heterogeneous and, in a mechanical model, the stress tensor is not any more spherical 
as in homogeneous bulks. The aim of this article is to show that a square-gradient 
functional taking into account the volume liquid free energy corrected with two 
surface liquid density functionals is a mean field approximation allowing to study 
structures of very thin liquid nanofilms near plane solid walls. The model determines 
analytically the concept of disjoining pressure for liquid films of thicknesses of a very 
few number of nanometers and yields a behavior in good agreement with the shapes 
of experimental curves carried out by Derjaguin and his successors. 
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1 Introduction 



The technical development of sciences allows us to observe phenomena at 
length scales of a very few number of nanometers. This nanomechanics infers 
applications in numerous fields, including medicine and biology. It reveals 
new behaviors, often surprising and essentially different from those that are 
usually observed at macroscopic and also at microscopic scales [1]. Currently 
simple models proposing realistic qualitative behaviors need to be developed 
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in different fields of nanoscicnccs even if their comparison with experimental 

data may be criticized at a quantitative level. 

As pointed out in experiments with water, the density of liquid water is found 
to be changed in narrow pores [2]. The first reliable evidence of this effect 
was reported by V.V. Karasev, B.V. Derjaguin and E.N. Efremova in 1962 
and found after by many others ([3], pages 240-244). In order to evaluate 
the structure of thin interlayers of water and other liquids, Green-Kelly and 
Derjaguin employed a method based on measuring changes in birefringence 
[4] ; they found significant anisotropy of water interlayers. 
In a recent article, the equations of motion of thin films were considered by 
taking into account the variation of the disjoining pressure along the layer 
[5]. The aim of this paper is to study, by means of a continuous mechanical 
model, the disjoining pressure and the behavior for very thin liquid films at 
the mesoscopic scale of a few number of nanometers. 

Since van der Waals at the end of the nineteenth century, the fiuid inhomo- 
geneities in liquid-vapor interfaces were represented in continuous models by 
taking into account a volume energy depending on space density derivative 
[6,7,8,9,10]. Nevertheless, the corresponding square-gradient functional is un- 
able to model repulsive force contributions and misses the dominant damped 
oscillatory packing structure of liquid interlayers near a substrate wall [11,12]. 
Furthermore, the decay lengths are correct only close to the liquid-vapor criti- 
cal point where the damped oscillatory structure is subdominant [13]. In mean 
field theory, weighted density-functional has been used to explicitly demon- 
strate the dominance of this structural contribution in van der Waals thin films 
and to take into account long-wavelength capillary- wave fluctuations as in pa- 
pers that renormalize the square-gradient functional to include capillary wave 
fluctuations [14]. In contrast, fluctuations strongly damp oscillatory structure 
and it is mainly for this reason that van der Waals' original prediction of a 
hyperbolic tangent is so close to simulations and experiments [15] . 
To propose an analytic expression in density-functional theory for liquid film 
of a very few nanometer thickness near a solid wall, we add a liquid density- 
functional at the solid surface and a surface density functional at the liquid- 
vapor interface to the square-gradient functional representing the volume free 
energy of the fiuid. This kind of functional is well-known in the literature 
[16]. It was used by Cahn in a phenomenological form, in a well-known pa- 
per studying wetting near a critical point [17]. An asymptotic expression is 
obtained in [18] with an approximation of hard sphere molecules and London 
potentials for liquid-liquid and solid-liquid interactions: in this way, we took 
into account the power-law behavior which is dominant in a thin liquid film 
in contact with a solid. 

The disjoining pressure 11 is a well adapted tool for a very thin liquid film of 
thickness h. In cases of Lifshitz analysis [19] and van der Waals theory, the 
disjoining pressure behaviors are respectively as 11 ~ h^^ and 11 ~ exp(— /i). 
None of them represents correctly experimental results for a film with a thick- 
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ness ranging over a few nanometers. 

Then, the gradient expansion missing the physically dominant damped oscilla- 
tory packing structure of the liquid near a substrate wall and only working up a 
smooth exponential decay is corrected by surface energies issued from London 
forces which model power-law dispersion interaction: since the only structure 
that the square gradient functional can yield is monotonic exponential, the 
surface energies take account of the eventual dominance of attractive power- 
law dispersion interactions. In fact power-law wings are physically present for 
liquid film of several nanometers and it is the reason we propose a study only 
for films in a range of few nanometers. 



2 The density-functional 



In our model, the free energy density-functional of an inhomogeneous liquid 
in a domain O of boundary dO is taken in the general form 

where e is the specific free energy and ip is a. generic surface free energy of dO. 
In our problem, we consider a horizontal plane liquid layer (L) contiguous 
to its vapor bulk and in contact with a plane solid wall (S); the z-axis is 
perpendicular to the solid surface [S). The liquid film thickness is denoted 
by h. Far from its critical point, the liquid at level z = h is situated at a 
distance order of two molecular diameters from the vapor bulk and the liquid- 
vapor interface is assimilated to a surface (E) at z — h. Then, the free energy 
density-functional (1) gets the particular form 

F ^ I I I e dv + [ [ (I) ds + [ [ ds. (2) 

J J J{L) J J{S) J JCE) 

where (p is shared in two parts 4> and ip respectively associated with {S) and 

• In Rel. (2), the first integral (energy of volume (L)) is associated with 
square-gradient approximation when we introduce a specific free energy of the 
fluid 

e = e{p, (3) 

at a given temperature ^ as a function of density p and f3 = (gradp)^. Spe- 
cific free energy e characterizes both fiuid properties of compressibility and 
molecular capillarity of liquid- vapor interfaces. In accordance with gas kinetic 
theory [20], scalar A = 2pe^jj{p,l3) (where e^jj denotes the partial derivative 



3 



with respect to /3) is assumed to be constant at a given temperature and 

p£ = Pq;(p) + ^ (grad p)^ (3) 

where the term (A/2) (grad pY is added to the volume free energy pa{p) of a 
compressible homogeneous fluid. We denote the pressure term associated with 
specific free energy a{p) by 

P(p)=A'(p). (4) 

• In Rel. (2), the second integral (energy of surface (S)) is defined through 
a model of molecular interactions between the fluid and the solid wall. In 
fact, near a solid wall, the London potentials of liquid-liquid and liquid-solid 
interactions are 

ipii — — ^ , when r > ui and — oo when r < ai , 
fis ~ ' when r > 5 and ipi^ — oo when r < 5 , 

where Cu and Qg are two positive constants, ai and (7^ denote liquid (fluid) 
and solid molecular diameters, 6 = \{cri + cTs) is the minimal distance be- 
tween centers of hquid and sohd molecules [21]. In the theory of additive 
and non-retarded molecular interactions, coefficients cu and are connected 
with Hamaker constants An and Aig through the relations An = n'^ciipf and 
Ais = tt'^CisPiPs, where pi and Pg respectively denote liquid bulk and solid 
densities [22]. Forces between liquid and solid have short range and can be 
simply described by adding a special energy at the surface. This is not the 
entire interfacial energy: another contribution comes from the distortions in 
the liquid density profile near the wall [18,23]. Finally, for a plane solid wall 
(at a molecular scale), this surface free energy is obtained in the form 

0(P) = -7iP+ ^72P^- (5) 

Here p denotes the liquid density value at the wall. The constants 7^, 72 are 

positive and given by the relations 7i = 2 Psj I2 — 2 — o > where 

125 mimg 125 

mi et rus respectively denote the masses of liquid (fluid) and solid molecules 
[18]. Moreover, we have A = ^ . 

3(7; mf 

• In Rel. (2), let us consider the third integral. The conditions in the vapor 
bulk are grad p = and Ap = with A denoting the Laplace operator. Far 
from the critical point, a way to compute the total free energy of the complete 
liquid-vapor layer is to add the energy of the liquid layer (L) located between 
2; = and z = h (first integral of Rel. (2)), the surface energy of the solid 
wall {S) at 2; = (second integral of Rel. (2)), the energy of the hquid- vapor 
interface of a few Angstrom thickness assimilated to a surface (S) at 2; = /i 
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and the energy of the vapor layer located between z = h and z = +00 [24]. 
The liquid at level z = h is situated at a distance order of two molecular 
diameters from the vapor bulk and the vapor has a negligible density with 
respect to the hquid density [25]. In our model, these two last energies can 
be expressed by writing a unique energy ip per unit surface located on the 
mathematical surface (E) at z — h : hj a calculation like in [18], we can write 
i/j in a form analogous to expression (5) and also expressed in [23] in the form 
il'{p) = — 75P + |74P^; but with a wall corresponding to a negligible density, 
75 ~ and the surface free energy ip is reduced to 

V'(P) = y P', (6) 

where p is the liquid density at level z = h and 74 is associated with a distance 
of the order of the fluid molecular diameter (when ai ~ S, then 74 ~ 72). 
Consequently, due to the small vapor density, the surface free energy ip is the 
same as the one of a liquid in contact with a vacuum and expressed by the 
third integral of Rel. (2). 

Such a form of density functional restricted to the first two integrals was pri- 
mary expressed by Cahn; Cahn's study used a graphic representation where 
energy integrals were presented as different areas in an energy-density plane 
[17]. Analytical computations were also tested in [24] but without taking ac- 
count of a complete volume free energy in form (3). 

With our previous functional approximation, we obtain the equations of equi- 
librium (or motion) and boundary conditions for a thin liquid film damping 
a solid wall. We can compute the liquid layer thickness. The normal stress 
vector acting on the wall remains constant through the layer and corresponds 
to the gas- vapor bulk pressure which is usually the atmospheric pressure. 
We obtain analytical results expressing the profile of density of very thin layer 
at a mesoscopic scale. We deduce an analytic expression of the disjoining pres- 
sure computed for different solid materials in contact with nanometer scale 
liquid layers. For all 1 know, such results have not been obtained in the litera- 
ture by using both a continuous mechanical model and a differential equation 
system. 

It is wondering to observe that the density-functional theory expressed by a 
simple model correcting van der Waals' one with surface density-functionals 
at the wall and the interface, enables to obtain a representation of the disjoin- 
ing pressure for very thin films which fits in with experiments by Derjaguin 
and others. This result is obtained without too complex weighted density- 
functionals and without taking account of quantum effects corresponding to 
an Angstrom length scale. So, this kind of functional may be a good tool 
to analytically study liquids in contact with solids at a very small nanoscale 
range. 
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3 Equation of motion and boundary conditions 



In case of equilibrium, functional F is minimal with respect to the vector fields 

of virtual displacement classically defined (as in [26]) and yields the equation of 
equilibrium of the inhomogeneous liquid and the boundary conditions between 
liquid, vapor and solid wall. In case of motions we simply add the inertial 
forces and the dissipative stresses in the equation of equilibrium (to refer to 
the well-known explicit calculations, see for example [5,27,28]). 



3. 1 Equation of motion 



The equation of motion is 

pT — div (<T + <T^) — p grad Q, 
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where T is the acceleration vector, Vt the body force potential, a the stress 
tensor generalization and cr^ the viscous stress tensor, 

cr — — pi — \ grad p ® grad p, 

where p = p^e p—p div (A grad p) is different from the pressure term P defined 
in (4). 

For a horizontal layer, in an orthogonal system of coordinates such that the 
third coordinate z is the vertical direction, all physical quantities in the layer 
depend only on z and the stress tensor cr of the thin film gets the form 
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Let us consider a thin film of liquid at equilibrium (gravity forces are ne- 
glected). The equation of equilibrium is 



div cr = . 

Equation (8) yields a constant value for the eigenvalue 03, 
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(8) 
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where P^^ denotes the pressure P^p^J in the vapor bulk, where p^^ is the den- 
sity of the vapor mother bulk bounding the liquid layer. Eigenvalues ai, 02 are 
not constant and depend on the distance z to the solid wall [3]. At equilibrium, 
general Eq. (7) is equivalent to 

grad(//(p)-AAp)=0, (9) 

where p, is the chemical potential at temperature 9 defined to an unknown 
additive constant [5,27]. The chemical potential is a function of P (and 9) but 

can be also expressed as a function of p (and 9). At temperature 9, we choose as 
reference chemical potential p^ = Pg{p), null for the bulks of densities pi and p^ 
of phase equilibrium. Due to Maxwell rule, the volume free energy associated 
with Pg is go{p) — Po where Po — P{pi) — P{Pv) is the bulk pressure and 
9o{p) — Ip Poip) is ^^^^ for the liquid and vapor bulks of phase equilibrium. 
The pressure P is 

P{p) = PPoip)-go{p) +Po- (10) 

Thanks to Eq. (9), we obtain in all the fluid and not only in the liquid layer, 

Poip) - ^^P^ PoiPb), 

where p^{pi^ is the chemical potential value of a liquid mother bulk of density 
such that PoiPb) = PoiPvb)- must emphasis that -P(Pb) and P{PyJ are 
unequal as for drop or bubble bulk pressures. The density pj, is not a fluid 
density in the layer but the density in the liquid bulk from which the layer 
can extend (this is why Derjaguin used the term mother liquid [3], page 32). 
In the liquid layer (L), 

^T^ = Pbip), with p6(p) = p„(p) - Po(p6). (11) 



3.2 Boundary conditions 



Condition at the solid wall (S) associated with the free surface energy (5) 
yields [28] 

where n is the external normal direction to the fluid. Equation (12) yields 

\ ' |z=0 

The sign of —j^ + 72 P|._„ determines the wettability of the fluid on the wall: 
the fluid damps the solid wall when 7^^ — ■j2P\^=o ^ ^ '^'^'^^ ^'^^ damp the 
solid wall when 7^ — ^2P\^=o ^ ^ [28,29]. 
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Condition at the liquid-vapor interface (E) associated with the free surface 
energy (6) yields 

Equation (14) defines a film thickness by introducing a reference point inside 
the liquid-vapor interface bordering the liquid layer with a convenient density 
at z = /i [24]. 

We notice that to study the stress tensor in the layer, we must also add to 
conditions (13,14) on density, the classical surface conditions on the stress 
vector associated with the total stress tensor <t + <t^ [28] . 



4 The disjoining pressure for horizontal liquid films 



We consider fluids and solids at a given temperature 6. The hydrostatic pres- 
sure in a thin liquid layer located between a solid wall and a vapor bulk differs 
from the pressure in the contiguous liquid phase. At equilibrium, the addi- 
tional pressure in the layer is called the disjoining pressure [3]. Clearly, the 




H 



Fig. 1. Diagram of the technique for determining the disjoining pressure isotherms 
of wetting films on a solid substrate: a circular wetting film is formed on a flat 
substrate to which a microporous filter with a cylindrical hole is clamped. A pipe 
connects the filter filled with the liquid to a reservoir containing the liquid mother 
bulk that can be moved by a micrometric device. The disjoining pressure is equal to 
n = (P6 - P„J gH (Prom Ref [3], page 332). 

disjoining pressure could be measured by applying an external pressure to 
keep the layer in equilibrium. The measure of the disjoining pressure is either 
the additional pressure on the surface or the drop in the pressure within the 
mother hulks that produce the layer. In both cases, the forces arising during 
the thinning of a film of uniform thickness h produce the disjoining pressure 
n(/i) of the layer with the surrounding phases; the disjoining pressure is equal 
to the difference between the pressure P^^ on the interfacial surface (which is 
the pressure of the vapor mother bulk of density p^J and the pressure in 
the liquid mother bulk (density p^) from which the layer extends: 
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Fig. 2. Experimental curves of the disjoining pressure (issued from Ref. [31]): the 
h-axis unit is the nanometer; the H-axis unit is about some atmospheres. The rectan- 
gle area corresponds to a domain where the curves of the disjoining pressure ( which 
depends on the quality of the wall) do not have a behavior as h~^ like in classical 
theory of thin films [19]. This experimental behavior corresponds to nanofilms of a 
very few number of nanometers. This behavior is similar to the one obtained from 
our model of functional which takes into account both wall and liquid-vapor interface 
effects. 



The most classical apparatus to measure the disjoining pressure is due to She- 
ludko [30] and is described on Fig. (1). The film is so thin that the gravity 
effect is neglected across the layer. Experimental curves of the disjoining pres- 
sure were first obtained by Derjaguin. The behavior of the disjoining pressure 
for a nanofilm in [31] seems strongly different from the one obtained for thin 
hquid film in [19] (see Fig. 2). 

If gb{p) = go{p) - go{Pb) - Po{Pb){P - Pb) denotes the primitive of /if,(p), null 
for p^, we get from Eq. (10) 

n(Pfe) = -9b{pvJ, (15) 
and an integration of Eq. (11) yields 

The reference chemical potential linearized near pi {respectively p^) is p^ip) = 
c2 / \ 

— ip — pi) respectively Poip) = ~iP " Pv) where q (respectively Cy) is the 
Pi \ Pv J 

isothermal sound velocity in liquid bulk pi {respectively vapor bulk p^) at 
temperature 9 [32] . In the liquid and vapor parts of the liquid- vapor film, Eq. 
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(11) yields 



= -j{p - Pb) (liquid) and A— = -^{p - p^J (vapor). 
The values of Po{p) are equal for the mother densities p^^ and p^, 

-z(Pb-Pi) = /^o(Pfe) = l^oiPvJ = —{Pv, - Pv): and consequently, 

Pv 



Pi 



Cf {Pb-PlY 



Pi 



In the liquid and vapor parts of the complete liquid- vapor layer we get the 
first expansion of the free energy, null when p — Pi and p — Py respectively, 

9o{p) = 7^{p - Pif (Hquid) and g„{p) ^ -^{p - p^f (vapor). 

^Pl ^Pv 

From definition of gh{p) and Eq. (15) we deduce the disjoining pressure 



^{P,) = ^{Pi-P,) 
^Pi 



Pi + Pft - P J 2 + 



c Upb-piY 

cl Pi . 



(17) 



/ (? { V 

Far from the critical point, due to I 2 + — ^ ^— ) <^ pi + p^, we get 



Pi 



n(p, 



'^Pi 



(Pi ~ Pb)- Now, we consider a film of thickness h; the density 



profile in the liquid part of the liquid- vapor film is solution of the differential 
equation, 

(18) 



dz"^ Pi 

. .dp ^ .dp 

with A— = -7i + 72 Pi^_Q and A— = -74 p, , 



With defining r such that t — cij \j\pi , where d = 1/r is a reference length 
and 73 = At, the solution of Eq. (18) is 



P = P6 + Pie ^^ + P2e^^ 



(19) 



where the boundary conditions at 2; = and z — h yield the values of p^ and 
P2 satisfying 

(72 + 73)Pi + (72 - 73)^2 = 7i - 72P6, 
_e-^^(73 - ^Jp^ + £^^^(73 + 74)p2 = -74P6 . 

The hquid density profile is a consequence of Eq. (19) when z e [0,h], 
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I (7l-72P&)(73 + 74)6^^ + (72 -73)74^6 ^-rz ^ 

^ (72 + 73) (73 + 74) e'*^ + (73 - 74) (72 - 73)6"''^ 



-(72 + 73)74Pb + (7i - 72P&)(73 - 74 je 
(72 + 73) (73 + 74)e''^ + (73 - 74) (72 - 73)e 



-hr 

(20) 



-hr 



Equations (16,19) together with gb{p) = {cf/2pi){p — p^)^ for the hquid part 
of the layer yield 

n(pJ = -^PiP2- (21) 

Pi 

The disjoining pressure is an invariant through the liquid film and its value is 
function of both and h, 



2cf 



^(Pb)-— (7i-72P6)(73 + 74)e"" + (72 
Pi ^ 



ish^Pb 



X 



(72 + 73)74Pb - (7i - 72P6)(73 - 74)e' 



-hr 



[(72 + 73) (73 + 74)e^^ + (73 - 74) (72 - 73)e ''l^ 



(22) 



By identification of expressions (17) and (22), we get a relation between h and 
Pjj. Consequently, we get a relation between the disjoining pressure Il{Pb) and 
the thickness h of the liquid film. For the sake of simplicity, wc denote the 
disjoining pressure as a function of h at temperature by 11 = Il{h). 
In experiments, for liquid in equilibrium with bubbles - even with a bubble 
diameter of a few number of nanometers - we have pj, ~ pi [33] . Consequently, 
the disjoining pressure is expressed as a function of h in the approximative 
form 



Pi 



r 

U{h) = -i [(7i - 72P;)(73 + 74)e'^" + (72 - 73)74^; 

(72 + 73)74^ - (7i - 72A)(73 - 74)e"''^ 



X 



[(72 + 73) (73 + 74)e''^ + (73 - 74) (72 - 73)e ''I 



Let us notice an important property of the mixture of a fluid far under its 
critical point and a perfect gas, where the total pressure is the sum of the 
partial pressures of the components [32]: at equilibrium, the partial pressure 
of the perfect gas is constant through the liquid-vapor-gas layer where the 
perfect gas is dissolved in the liquid. The disjoining pressure of the mixture 
is the same as for a single fluid and calculations and results are identical to 
those previously obtained. 
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Physical constants 


Cll 


<7i 


nil 


Pi 


ci 


Water 


1.4 X 10-5^ 


2.8 X 10-s 


2.99 X 10-23 


0.998 


1.478 X 10^ 


Physical constants 


Cls 




rris 




(5 


Silicon 


1.4 X 10-5^ 


2.7 X 10"^ 


4.65 X 10-23 


2.33 


2.75 X 10-^ 


Deduced constants 


A 


7i 


72 = 74 


73 


d 


Results (water-silicon) 


1.17 X 10"^ 


81.2 


54.2 


506 


2.31 X 10-^ 



Table 1 



The physical values associated with water and silicon are obtained in references 
[21,34] and expressed in c.g.s. units (centimeter, gramme, second). No information 
is available for water-silicon interactions; we assume that cn = cig. The deduced con- 
stants are obtained from physical values by means of formulae obtained in sections 
2-4. 

5 A comparison of the model with experiments 

Our aim is not to propose an exhaustive study of the disjoining pressure 
for all physicochemical conditions associated with different fluids bounded 
by different walls, but to point out an example such that previous modeling 
appropriately fits with experimental data. At ^ = 20° Celsius, we successively 
consider water wetting walls (a wall of silicon is the reference material) and 
water not wetting a wall. 

Due to Eq. (20), Fig. 3 represents water liquid density profiles in the nanolayer. 
We verify the consistency of the model: 

• The density gradient is large at a few nanometer range from the solid wall 
and consequently in this domain, the liquid is inhomogeneous, 

• The boundary condition (14) is well adapted to our model of functional 
and determines the position where the phase transition between liquid and 
vapor occurs: condition (14) yields a density value of the fluid corresponding 
to an intermediate density which can be associated with a dividing surface 
separating liquid and vapor in the liquid-vapor interface. Due to the film 
instability, we will see further down that graph (b) in Fig. 3 is unphysical. 

We have drawn disjoining pressure profiles deduced from analytical expressions 
given in section 4; the graphs relate to Rel. (22). 

Graphs are associated with several cases when water damps the solid wall (at 
the wall, the water density is closely p; and 7^ — 72P; > 0) and a case when 
water does not damp the solid wall (7^ — ^2Pi < 0)- 

According to different physical values, graphs (a), (b) and (c) of Fig. (4) 
represent the disjoining pressure profiles for water in contact with a plane 
solid wall at 20° Celsius. 
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Fig. 3. Graphs of liquid density profiles in a nanofilm. Results are given in two cases: 

Graph (a) corresponds to the case of a liquid damping a solid wall ('Ji — < 
wetting data are associated with liquid water on silicon (cf. graph (h) in Fig. 4)- 
Graph (b) corresponds to the unphysical case of non-wetting liquid ('Ji — ^2Pi > 
non-wetting data are associated with liquid water and a wall such that 7^ = 30, all 
the other deduced constants A, 72, 73, 74 being unchanged (cf. graph (d) in Fig. 4)- 
The unit of the x—axis is 5 = 2.75 x 10^^ cm (2.75 Angstrom), the unit of the 
y—axis is the liquid water density at 2(f Celsius (approximatively 1 g/crr?). 



Values of 71 


110 


81.2 


58 


30 


Corresponding graphs in Fig. 4 


(a) 


(b) 


(c) 


(d) 



Table 2 

The numerical data A, 72, 73, 74 corresponding to liquid water are unchanged. They 
have the values presented in Table 1. Only the values of 7^^ (in c.g.s. units) depending 
on the behavior of the solid walls are different following graphs (a), (b), (c), (d) 
drawn in Fig. 4. 



Graph (b) corresponds to a silicon solid wall. Graphs (a) and (c) respectively 
correspond to water wetting more strongly the wall than a silicon wall and 
water wetting less strongly the wall than a silicon wall. 

In [5], we studied the stability of nanofilm. In accordance with results in [3], 
Graph (a) is associated with a stable nanolayer for any liquid film thickness 
because for all h, dU.{h, 9)/dh < 0. In graphs (b) and (c), values of h for which 
the liquid nanolayer is stable correspond to a domain where (9n(/i, 9)/dh < 0, 
corresponding to h values greater than a particular value hs depending on 7^. 
Graph (d) differs from previous ones as that liquid water does not damp the 
solid wall. The graph corresponds to an unstable nanolayer and does not exist 
physically. In the non-wetting case, liquid nanolayers are unstable and they 
are associated with compression instead of suction in experiments by Sheludko 
[30]. 

We notice that graphs (a), (b), (c) in Fig. 4, experimental graph in Fig. 2 and 
graphs in experimental literature (as in [3,31]) exhibit quite similar behaviors. 
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Fig. 4. Several graphs represent the disjoining pressure for liquid water at 20° C 
in contact with different plane solid walls: graphs (a), (h), (c) correspond to water 
damping the wall; graph (d) corresponds to an unphysical case of water not damping 
the wall. The unit of the x—axis is S, the unit of the y—axis is 10^ Barye ( one 
atmosphere) . 



6 Conclusion 



We have studied liquid nanofilms in contact with plane solid walls. For layer 
thicknesses of some nanometers, the theoretical graphs of the disjoining pres- 
sure correctly draw the behavior of experiments by Derjaguin and others [3,30]. 
The proposed analytical method is different from Lifschitz one in which layers 
were considered with uniform density liquids [19]. In our model correspond- 
ing to thin liquid nanofilms, liquids are considered as inhomogeneous near 
the solid walls. The density distribution in liquid nanofilms depends on the 
physico chemical characteristics of walls: when the liquid damps the wall, we 
have an excess of fluid density at the wall and the fluid is denser at the wall 
than in the liquid bulk; the contrary happens when the liquid does not damp 
the wall. 

These analytical results and the liquid density profiles are in accordance with 
experimental works by Derjaguin and others [3,4,11]. 
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